OncoTargets and Therapy 



Q Open 



Access Full Text Article 



Dovepress 

open access to scientific and medical research 

ORIGINAL RESEARCH 



Expression profiles and initial confirmation 
of long noncoding RNAs in Chinese patients 
with pulmonary adenocarcinoma 



Xin Zhao* 
Wen Zhu* 
Wangjian Zha 
Feifei Chen 
Zhenzhen Wu 
Yanan Liu 
Mao Huang 

Department of Respiratory Medicine, 
The First Affiliated Hospital, Nanjing 
Medical University, Nanjing, People's 
Republic of China 

*These authors contributed equally to 
this study and share first authorship 



Correspondence: Mao Huang 
Department of Respiratory Medicine, The 
First Affiliated Hospital, Nanjing Medical 
University, 300 Guangzhou Road, Nanjing 
210029, People's Republic of China 
Tel +86 25 6813 6269 
Fax +86 25 8367 3567 
Email hm6l I4@l26.com 



This article was published in the following Dove Press journal: 
OncoTargets and Therapy 
I July 2014 

Number of times this article has been viewed 



Background: The purpose of this study was to investigate differentially expressed long 
noncoding RNAs (IncRNAs) in pulmonary adenocarcinoma tissue and adjacent noncancerous 
tissue from Chinese patients using IncRNA expression microarray and preliminary analysis. 
Methods: RNA extracted from three paired pulmonary adenocarcinoma tissue and adjacent 
noncancerous tissue specimens was used to synthesize double-stranded complementary DNA 
after labeling and hybridization. The complementary DNA was labeled and hybridized to 
the IncRNA expression microarray, and array data were analyzed for hierarchical clustering. 
Gene coexpression networks were constructed to identify interactions among genes. To validate 
the microarray findings, we measured the relative expression levels of four random differen- 
tially expressed IncRNAs in the same tissue used for microarray using real-time quantitative 
polymerase chain reaction. The expression level of one IncRNA, AK1 24939, in the paired 
pulmonary adenocarcinoma/adjacent noncancerous tissue of another 30 patients was measured 
using real-time quantitative polymerase chain reaction. The experimental data were further 
analyzed and compared with clinical features. 

Results: Of 39,000 IncRNAs investigated, 704 were differentially expressed in pulmonary 
adenocarcinoma tissue; 385 were upregulated and 319 were downregulated compared with 
those in the adjacent noncancerous tissue (fold change ^2 and £- 2, P<0.05). AK124939 
expression levels in poorly differentiated adenocarcinoma tissue were lower than those found 
in well to moderately differentiated adenocarcinoma tissue (P=0.05). 

Conclusion: There are significant differences in the IncRNA expression profiles in Chinese 
patients with pulmonary adenocarcinoma. LncRNAs such as AK124939 may be anticancer 
factors related to the progression of pulmonary adenocarcinoma. 
Keywords: pulmonary adenocarcinoma, long noncoding RNA, microarray 

Introduction 

One of the greatest surprises in modern biology was the discovery that the human 
genome encodes only about 20,000 protein-coding genes, representing less than 2% 
of the total genome sequence. 1 However, the advent of tiling resolution genomic 
microarrays and whole-genome and transcriptome sequencing technologies revealed 
that at least 90% of the genome is actively transcribed. 2,3 It was found that the 
complexity of the human transcriptome exceeds that of a collection of protein-coding 
genes and their splice variants; there is extensive antisense, overlapping, and noncoding 
RNA expression. 4-6 

Although initially considered spurious transcriptional noise, recent evidence 
suggests that the so-called genomic "dark matter" may play a major biological role 
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in cellular development and metabolism. 7 " 10 The newly 
discovered long noncoding RNA (IncRNA) genes are one 
such element, with developmental and tissue-specific expres- 
sion patterns and aberrant regulation in a variety of diseases, 
including cancer. 11 " 19 

Discovery of IncRNAs 

LncRNAs were first described in a rat full-length comple- 
mentary DNA library. LncRNAs were originally isolated in 
proximity to the noncoding RNA of protein-coding genes 
and were classified as sense, antisense, bidirectional, intronic, 
and intergenic according to their positions relative to the 
protein-coding genes in the genome. Initially, IncRNAs were 
considered transcriptions of genome "noise", a byproduct 
of RNA polymerase II low-fidelity transcription with no 
biological function, 9 and were thus neglected. In recent 
years, gradual recognition of the role of IncRNAs in X 
chromosome silencing, genomic imprinting, chromatin 
modification, transcription, transcription interference, and 
nuclear transport regulation have led to IncRNAs becoming 
a focus of research. 

Biological characteristics and function of 
IncRNAs 

LncRNAs are transcriptions of >200-nucleotide RNAs and 
are typically transcribed by RNA polymerase II. Besides 
temporal and spatial expression specificity, they lack mean- 
ingful open reading frames and do not code for proteins, 
instead directly regulating gene expression in the form of 
RNA at various levels, and are related to the occurrence and 
development of many diseases, including cancer. 8 

LncRNA genes are widely distributed in the genome and 
might be located in messenger RNA (mRNA) gene exons 
or introns, or between mRNA genes. At the same time, sub- 
cellular localization of IncRNAs is complex, as they can be 
located in the nucleus or cytoplasm. 10 LncRNAs and mRNAs 
share many similar characteristics and overlapping areas: the 
genome "transcription hotspot" usually can be transcribed 
into both mRNA and IncRNA; each of the two chains of 
protein-coding genes can be transcribed into IncRNA; 2 the 
same protein-coding gene can be transcribed into different 
mRNAs or IncRNAs; and a part of IncRNAs and mRNAs 
can be formed by alternative exon splicing combinations 
with polyadenylic acid in the 3' terminal. 

Generally, IncRNAs are mutable, but many contain highly 
conserved elements, which may be the result of a rapid 
adaptation to evolution 20 that differs from that of mRNAs, 
which must ensure codon correctness to prevent open reading 



frame change. LncRNAs can only maintain a high level of 
conservation in sections that sustain secondary structure sta- 
bility or those playing a key role in the functions of IncRNA. 
Studies have shown that mammalian IncRNAs are poorly 
conserved at the primary structure level while being highly 
conserved at secondary structure level. 6,21 

How IncRNAs regulate gene expression varies widely 
and is manifested in diverse IncRNA regulatory mechanisms. 
The effect of IncRNAs varies with the mechanism involved. 
LncRNAs can function as a main transcription factor, or 
a common regulatory factor that plays a role with other 
components. Gene expression regulation can generally be 
divided into: 

1 . Regulation at the epigenetic modification level through 
interaction with various chromatin-modifying enzymes 
that modify chromatin by changing its conformation or 
activity, or repressing expression of the genes responsible. 
At the embryonic stage of development in particular, 
IncRNAs are involved in late silencing of allelic gene 
expression, maintaining epigenetic status, which is 
essential for normal development and cell differentiation 
in multicellular animals. 

2 . Regulation at the transcriptional level through regulation 
of the combination and assembly of transcription factors, 
combination with regulatory DNA sequences to form 
three chain complexes, regulation of RNA polymerase 
II transcription, and transcription interference. 

3 . Regulation at the post-transcriptional level through for- 
mation of double-stranded RNA with complementary 
mRNA, affecting mRNA processing, splicing, transporta- 
tion, translation, and degradation. 

Emerging role of IncRNAs 
in tumorigenesis 

It has long been suspected that differential IncRNA expres- 
sion plays a role in cancer; however, this suspicion has lacked 
strong supporting evidence. With advances in cancer tran- 
scriptome profiling and accumulating evidence supporting 
the function of IncRNAs, a number of differentially expressed 
IncRNAs have been associated with cancer. LncRNAs have 
been implicated in the regulation of a range of biological 
functions and the disruption of some of these functions, such 
as genomic imprinting and transcriptional regulation, thus 
playing a critical role in development of cancer. 

Lung cancer is a malignant tumor with the highest morbid- 
ity and mortality worldwide. In 2003, Ji et al used sequenc- 
ing and quantitative reverse transcription-polymerase chain 
reaction (qRT-PCR) and found that among 225 non-small 
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cell lung cancer (NSCLC) patients, metastasis-associated 
lung adenocarcinoma transcript 1 (MALAT 1 ) IncRNA was 
overexpressed in 70 samples from patients with metastatic 
tumors. Further, the overexpression was stage-specific and 
tissue-specific, suggesting a close relationship between 
MALAT 1 and the high metastatic potential and poor prog- 
nosis of NSCLC, thus implicating MALAT1 as a prognostic 
indicator of survival in patients with stage I NSCLC. 22 The 
MALAT1 transcript is a >8,000-nucleotide noncoding RNA 
located on chromosome llql3. Although it is expressed 
in normal human tissue and various tumors such as breast 
cancer, bladder cancer, liver cancer, and pancreatic cancer, 
overexpression in NSCLC is the most significant. Using 
cell migration experiments and qRT-PCR, Tano et al found 
that MALAT 1 -silencing mediated by small interfering RNA 
regulated the transcription and post-transcriptional expres- 
sion of migration-related genes to affect cell migration ability 
in pulmonary adenocarcinoma, disrupting the migration of 
adenocarcinoma cells in vitro. 23 These studies indicate that 
MALAT 1 regulates metastasis. In addition, using overexpres- 
sion and small interfering RNA interference methods, Trip- 
athi et al found that MALAT 1 and serine/arginine-rich (SR) 
proteins interact to affect the distribution of SR and other 
splicing factors in small nuclear areas. MALAT 1 affected 
alternative splicing of the pre-mRNA precursor by regulating 
the level of phosphorylated SR protein, 24 which is of great 
significance in the study of SR protein splicing regulation, 
although its role in tumorigenesis remains unknown. 

An important lung cancer, pulmonary adenocarcinoma is 
one of the most common human cancers, with poor 5 -year 
survival rates worldwide. In this study, we profiled IncRNA 
expression in Chinese patients with pulmonary adenocarci- 
noma and performed a preliminary analysis of the experi- 
mental and clinical data. 

Materials and methods 

Samples 

Three paired pulmonary adenocarcinoma tissue (T 1 , T2, T5) 
and adjacent noncancerous tissue (Nl, N2, N5) specimens 
from patients with grade II tumors and without a history of 
chemotherapy and radiotherapy were used for the microar- 
ray study. The tumor (T) and normal (N) tissues were stored 
in liquid nitrogen. The fourth T/N tissue pair used in the 
microarray was a mixture of the first three pairs: T6 was a 
combination of Tl, T2, and T5, and N6 was a combination 
ofNl,N2, andN5. 

Another 30 paired pulmonary adenocarcinoma/adjacent 
noncancerous tissue specimens from patients who had 



undergone surgical therapy at the First Affiliated Hospital of 
Nanjing Medical University between July 20 1 2 and February 
2013 and for which the postoperative pathological diagno- 
sis was confirmed to be pulmonary adenocarcinoma were 
obtained for the qRT-PCR study. 25 26 All specimens were 
collected after patients had signed informed consent forms 
and following review by the hospital ethics committee. 

Microarray and computational analysis 

RNA extracted from the three paired tissue specimens 
was used to synthesize double-stranded complementary 
DNA after labeling and hybridization. The complementary 
DNA was labeled and hybridized to an Agilent Human 
IncRNA + mRNA array version 2.0 (Agilent Technologies, 
Santa Barbara, CA, USA) designed with four identical arrays 
per slide (4 x 180,000 format). Each array contained probes 
interrogating 39,000 human IncRNAs and 32,000 human 
mRNAs. Each RNA was detected twice by the probes. The 
array also contained 4,974 Agilent control probes. 

The array data were analyzed for summarization, normal- 
ization, and quality control using GeneSpring software ver- 
sion 1 1 .5 (Agilent). To select differentially expressed genes, 
we used threshold values of a ^2- and s-2-fold change and 
a Benjamini-Hochberg corrected P-value of 0.05. The data 
were log 2 transformed and median-centered by genes using 
the Adjust Data function of CLUSTER version 3.0 software 
(Michiel de Hoon, Human Genome Center, University of 
Tokyo, Tokyo, Japan), then further analyzed with hierarchi- 
cal clustering with average linkage. 27 Lastly, we performed 
tree visualization using Java Tree View (Stanford University 
School of Medicine, Stanford, CA, USA). 

The microarray data discussed in this paper have 
been deposited in the National Center for Biotechnology 
Information Gene Expression Omnibus and can be accessed 
through the series accession number GSE56850 ( http://www. 
ncbi, nlm.nih.gov/geo/query/acc. cgi?acc=GSE56850 ). 

Coexpression network 

Gene coexpression networks are used to identify interactions 
among genes 28 according to the normalized signal intensity of 
specific expressed genes. We constructed network adjacency 
between two genes, i and j, defined as a power of the Pearson 
correlation between the corresponding gene expression pro- 
files xi and xj. We then obtained a gene adjacency matrix, 
M (i,j), 29 by computing the profiles, visualized as a graph, the 
topological properties of which were examined. To construct 
a visual representation, only the strongest correlations (>0.99) 
were drawn in the renderings. Each gene corresponded to a 
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node in the gene coexpression networks. Two genes connected 
by a line indicated a strong correlation (ie, red lines denote 
positive correlation, blue lines denote negative correlation.). 
In network analysis, a degree is the simplest, most important 
measure of the centrality of a gene within a network and 
determines the relative importance. A degree is denned as 
the number of directly linked neighbors. 30 

qRT-PCR 

The qRT-PCR was performed using PrimeScript RT Master 
Mix (Takara, Otsu, Japan) and SYBR Green Mix (ABI, 
Foster City, CA, USA). Total RNA was extracted from tissues 
using TRIzol (Invitrogen, Carlsbad CA, USA) and purified 
with a mirVana miRNA isolation kit (Ambion, Austin, TX, 
USA) according to the manufacturer's protocol. RNA purity 
and concentration were determined from optical density 
OD260/280 readings using a spectrophotometer (NanoDrop 
ND-1000; Thermo Fisher Scientific, Wilmington, DE, USA); 
RNA integrity was determined by 1 % formaldehyde denatur- 
ing gel electrophoresis. The RNA was reverse-transcribed into 
complementary DNA by combining 4 uL 5x PrimeScript RT 
Master Mix, 2 [iL template RNA, and RNase-free distilled 
water to a total volume of 20 |iL in a microtube, heating at 
37°C for 15 minutes, then at 85°C for 5 seconds, followed 
by immediate cooling on ice. The purity and concentration of 
the complementary DNA were determined from OD260/280 
readings using the NanoDrop ND-1000 spectrophotometer. 
Each 10 uL reaction contained 5 uE 2x SYBR Select Master 
Mix, 0.2 |iL each forward and reverse primers, a complemen- 
tary DNA template (up to 1 00 ng of complementary DNA was 
used in each reaction for optimal performance), and 4.6 uL 
of RNase-free distilled water, and was transferred to wells in 
an optical plate. The plate was sealed with an optical adhe- 
sive cover and centrifuged briefly to spin down the contents 
and eliminate air bubbles before being placed in a StepOne 
thermocycler (Life Technologies, Carlsbad CA, USA). The 
thermal cycling conditions were as follows: uracil-DNA gly- 
cosylase activation at 50°C for 2 minutes, AmpliTaq DNA 
polymerase (Life Technologies) and primer activation at 
95°C for 2 minutes, denaturation at 95°C for 15 seconds, and 
annealing/extension at 60°C for one minute for 40 cycles. We 
programmed the instrument to perform a default dissociation 
step. We used glyceraldehyde-3-phosphate dehydrogenase as 
the internal control throughout the study. 

Statistical analysis 

All statistical analyses were performed using Statistical 
Package for the Social Sciences version 20 software 



(SPSS Inc., Chicago, IL, USA). We used one-way analysis 
of variance, Fisher's exact test, Pearson correlation analysis, 
and the two-tailed /-test (including the independent-samples 
f-test and paired /-test). Statistical significance was set at a 
P-valueof <0.05. 

Results 

LncRNA expression profiles in pulmonary 
adenocarcinoma 

Of 39,000 IncRNAs, 704 were differentially expressed in 
pulmonary adenocarcinoma tissue; 385 were upregulated 
and 319 were downregulated as compared with the adjacent 
noncancerous tissue (fold change s2 and ^-2, P<0.05). 
Hierarchical clustering revealed systematic variations 
between the expression of IncRNAs and protein-coding 
RNAs in the paired pulmonary adenocarcinoma/nontumor 
lung tissue specimens (Figure 1 A and B). The data indicate 
that a series of IncRNAs is frequently differentially expressed 
in pulmonary adenocarcinoma tissues. 

Bioinformatics analysis and validation 
of IncRNA expression profiles 

Gene coexpression networks are used to cluster thousands 
of transcripts into phenotypically relevant coexpression 
modules. 28 - 31 Given that coexpression networks may cor- 
respond to biological pathways 32 and many functions of 
protein-coding RNAs are present in the National Center 
for Biotechnology Information Reference Sequence data- 
base, we concentrated on coexpression with a high rate 
of protein-coding RNAs in the carcinoma coexpression 
network, selecting four IncRNAs (Table 1), ie, AK124939, 
ASO3690, ENST00000491282, and XLOC_000575, in 
pulmonary adenocarcinoma for characterization using this 
method. In the cancer coexpression network, AK124939 was 
connected to a protein-coding gene, PDZD2 (also known 
as AIPC, PIN1, PAPIN, PDZK3), which is involved in 
tumorigenesis and tumor development. ASO3690 and 
XLOC_000575 were each connected to two protein-coding 
genes; ENST00000491282 was connected to 17 protein- 
coding genes (Figure 2). 

To validate the microarray findings, we analyzed the 
expression of the abovementioned IncRNAs with qRT- 
PCR in the same three paired pulmonary adenocarcinoma 
tissue/adjacent noncancerous tissue specimens used for the 
microarray. Table 2 lists the corresponding primer sequences. 
The data confirmed that ASO3690, ENST0000049 1282, and 
XLOC_000575 were overexpressed in pulmonary adeno- 
carcinoma, whereas AK1 24939 expression was decreased 
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Figure I Differential expression of long noncoding RNAs in pulmonary adenocarcinoma tissue (Tl, T2, T5, T6) and adjacent noncancerous tissue (Nl, N2, N5, N6). 
Hierarchical clustering analysis of (A) 704 long noncoding RNAs and (B) 1 , 1 56 mRNAs that were differentially expressed (fold change >2, P<0.05). Red and green shades 
denote expression levels above and below the median value across all samples (-3 to +3), respectively. 
Abbreviations: IncRNA, long noncoding RNA; mRNA, messenger RNA; n, normal; t, tumor. 



(P<0.05 for all, Figure 3). These results were sufficiently 
consistent with that of the IncRNA expression profiles. 

AK 1 24939 may be related to progress 
of pulmonary adenocarcinoma 

To explore the correlation between IncRNA and pulmo- 
nary adenocarcinoma, we increased the number of samples 
studied by obtaining 30 paired pulmonary adenocarcinoma/ 
adjacent noncancerous tissue specimens and detected 
AK124939 and its related PDZD2 mRNA (Table 2 lists the 
primer sequences) using qRT-PCR. The transcript levels 
of AK 12493 9 and PDZD2 mRNA were both lower in pul- 
monary adenocarcinoma tissues (P<0.05) compared with 
the paired noncancerous lung tissues (Figure 4A and B). 
Pearson correlation analysis estimated ^0.741 (P<0.05) 
in the tumor tissues and ^0.905 (P<0.05) in the adjacent 
noncancerous tissues (Figure 5A and B). Therefore, there was 
a clear positive correlation between the IncRNA AK1 24939 
and PDZD2 mRNA. 

We then explored the relationship between level of 
AK 12493 9 expression and the clinical and pathological 
characteristics of patients with pulmonary adenocarcinoma. 
There was no significant correlation between AK124939 
expression and age, sex, tumor-node -metastasis (TNM) 
stage, tumor size, Karnofsky performance status, smoking 
history, or chemotherapy history. Although not statistically 



significant (P=0.05), the level of AK124939 expression was 
related to tumor differentiation grade. AK124939 expression 
in the poorly differentiated adenocarcinoma tissue was lower 
than that in well to moderately differentiated adenocarcinoma 
tissue (Table 3). 

Discussion 

Once considered transcription "noise", the possible role of 
IncRNAs in cancer development and occurrence has recently 
been revealed, for which studies are ongoing. Previous stud- 
ies reported more significant MALAT 1 IncRNA overexpres- 
sion in NSCLC compared with that in other cancers 22 and that 
IncRNA is a possible regulatory factor of metastasis. 23 Given 
these previous findings, we aimed to determine whether the 
same would hold true for other IncRNAs in pulmonary adeno- 
carcinoma, which has been described as the most prevalent 
form of NSCLC. 33 

We used a commercial microarray and qRT-PCR to 
investigate the differential expression of IncRNAs in pul- 
monary adenocarcinoma tissue specimens and compared 
them with that in adjacent noncancerous tissue specimens. 
Analysis of the microarray data revealed significant dif- 
ferential expression for 704 IncRNAs (385 upregulated and 
319 downregulated) in pulmonary adenocarcinoma tissue 
compared with that in adjacent normal tissue. Subsequent 
gene coexpression network analysis revealed that the 
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Figure 2 Long noncoding RNA subnetwork in the carcinoma coexpression network consisting of ASO3690, XLOC_000575, AK 1 24939, and ENST0000049 1 282 
(long noncoding RNAs, highlighted in yellow), and their direct neighbors. Genes highlighted in green are protein-coding RNAs. Red lines denote positive correlation; blue 
lines denote negative correlation. Node size represents the number of its neighbors. 



Table I The four studied long noncoding RNAs 



Database ID 


Accession 


P-value P-value Regulation 


Lnc 


Biotype-Transcripts_ Biotype-Gene_ 






corrected 


length 


Ensemble 68 Ensemble 68 


ASO3690 


None 


O.I 1 8582 0.000433 Up 


2,582 


Antisense Antisense 


XLOC_000S75 


None 


0.1 86826 0.0 1 9246 Up 


1,367 


LincRNA LincRNA 


AK 1 24939 


AK 1 24939 


0.1 44 1 22 0.00307 Down 


1,839 


Antisense Antisense 


ENST0000049I282 


None 


0.1 95805 0.02793 Up 


883 


Retained_intron LincRNA 


Abbreviations: Lnc, long noncoding; LincRNA, long intergenic noncoding RNA. 






Table 2 Primer sequences 








RNA 




Forward primer sequences 5'-3' 




Reverse primer sequences 5'-3' 


ASO3690 




CAAG C AAG AG AG CTGAGGATGA 




CTGGAGTCTCCTGTCTTCCCT 


AK 1 24939 




TGTGCATTGCCTCAGGTTCA 




AGCTGACTGGGAGTTGGTTC 


ENST0000049I282 




CGAAATGGATTCACGGTGCC 




TGCCAGATCAGGGTGTTGTC 


XLOC_000S75 




TGCCTTTGACCACATCCATCT 




ACCCAGTTCTAGCAGCCCAG 


PDZD2 




CCTGAACTTGGAAACGGCCAT 




AGTCACTCGATTGGAAGTCCG 


GAPDH 




CCACTCCTCCACCTTTGAC 




ACCCTGTTGCTGTAGCCA 


Abbreviations: LncRNA, long noncoding RNA; GAPDH, glyceraldehyde-3-phosphate dehydrogenase. 
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Figure 3 Quantitative reverse transcription-polymerase chain reaction validation 
of ASO3690, XLOC_000575, AK 1 24939, and ENST0000049I282 expression in the 
three paired pulmonary adenocarcinoma/nontumor lung tissue specimens used for 
the microarray. 



coexpression patterns of the IncRNAs and protein-coding 
RNAs of the pulmonary adenocarcinoma tissue differed from 
that in normal tissue. 

We characterized the IncRNAs, AK124939, ASO3690, 
ENST00000491282, and XLOC_000575, which had a high 
rate of protein-coding RNAs in the carcinoma coexpression 
network, and found that all but AK124939 were upregulated 
in pulmonary adenocarcinoma tissue. Similarly, mRNA 
expression of the tumor suppressor gene PDZD2, related 
to AK124939, was downregulated, indicating a correlation 
between the two. It has been demonstrated that PDZD2 is 
upregulated in primary prostate tumors and may be involved 
in the early stages of prostate tumorigenesis. 34,35 In light of 
these findings, more research is warranted to explore the 



difference in the mechanisms and relationship between 
IncRNAs and protein-coding RNAs. 

Every pattern of pulmonary adenocarcinoma is assigned a 
grade related to its metastatic or recurrence potential: grade 1 , 
lepidic; grade 2, acinar or papillary; and grade 3, solid or 
micropapillary, which correspond to low, intermediate, and 
high metastatic potential, respectively. Well differentiated 
tumors (histologic score 3, lepidic pattern) are associated with 
low rates of recurrence, while poorly differentiated tumors 
(histologic score of 5 and 6, solid or micropapillary pattern) 
are associated with higher rates of recurrence. Moderately 
differentiated tumors with a histologic score of 4 have pure 
acinar, papillary, or combined acinar and papillary as the pre- 
dominant and secondary predominant patterns. 36 - 37 Among the 
30 cases, four (13.3%) were well differentiated (grade 1/G1), 
ten (33.3%) were moderately differentiated (grade 2/G2), and 
16 (53.3%) were poorly differentiated (grade 3/G3). 

Further investigation of AK1 24939 expression in the 
30 paired pulmonary adenocarcinoma/adjacent noncancer- 
ous tissues revealed that its expression levels in the poorly 
differentiated adenocarcinoma tissue were lower than that 
in well to moderately differentiated adenocarcinoma tissue 
(albeit P=0.05). Given the low number of cases and lack 
of proportion, we merged Gl and G2 into a well to mod- 
erately differentiated group and compared these with the 
poorly differentiated (G3) group; there were no significant 
differences between Gl or G2 and G3. It is known that the 
poorer the tumor cell differentiation, the higher the degree 
of malignancy. In patients with the same TNM stage, his- 
topathological grade tends to have an effect on prognosis. 
It has been shown that IncRNAs are significant regulators 
of gene expression that interact with the major pathways of 
cell growth, proliferation, differentiation, and survival. 38,39 



B 



c3l°n 
'</> 

s 

Q. 

X 
CD 
CD 
> 

ra 

CD 



o 



P=0.049 



P=0.002 



o<5 



.9" 



| : : : : : :1 Nontumor 
Tumor 



c 15- 

o 

'in 

m 

P 

Q. 10 

x 

CD 
CD 
> 

1 5 

CD 
O) 




Nontumor 
Tumor 




Figure 4 Quantitative reverse transcription-polymerase chain reaction analysis of AK 1 24939 and PDZD2 expression in 30 paired pulmonary adenocarcinoma/adjacent 
noncancerous lung tissue specimens. (A) Paired t-test revealing statistically significant AK 1 24939 and PDZD2 downregulation in the 30 paired samples (P<0.05). (B) Lines 
connecting the level of AK 1 24939 and PDZD2 expression in each paired sample from the same patients. 
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Figure 5 Pearson correlation analysis estimated that there was a positive correlation between the IncRNA AK 1 24939 and PDZD2 mRNA. 
Notes: (A) f=0.74 1 (P<0.05) in the tumor tissues and (B) (=0.905 (P<0.05) in the adjacent noncancerous tissues. 
Abbreviations: IncRNA, long noncoding RNA; mRNA, messenger RNA. 



Alterations in the function of IncRNAs promote tumor 
formation and progression, and metastasis; for example, 
noncoding Nras functional RNA is involved in malignant 
transformation of colonic epithelial cells. 40 Therefore, if 
future studies involving a larger sample size confirm that 
there is a significant difference (P<0.05), we may assume 
that AK124939 is related to cell differentiation and tumor 
formation as a suppressor gene. From this perspective, 
multicenter, randomized concurrent controlled studies are 
required to determine the potential of AK 12493 9 as a diag- 
nostic and prognostic marker. Our findings indicate that, 
like MALAT1 IncRNA, AK124939 plays a potential role in 
targeted treatment of lung cancer, or as a marker of diagnosis 
or prognosis, and further research is required. 



Table 3 Patient characteristics 



Clinical characteristics 


n 


P-value 
for tumor 
tissues 


P-value for 
nontumor 
tissues 


Age ==60/>60 years 


14/16 


0.065 


0.100 


Male/female 


19/1 1 


0.248 


0.367 


Smoking history/ 


20/10 


0.124 


0.896 


no smoking history 








TNM l-ll/TNM III — IV 


21/9 


0.972 


0.487 


Tl (tumor size <3 cm)/T2 


13/17 


0.710 


0.537 


or T3 (tumor size >3 cm) 








NO (no LN metastasis)/N 1 


IS/IS 


0.474 


0.754 


or N2 (LN metastasis) 








KPS 80 or 90/ 100 


21/9 


0.826 


0.373 


Well/moderately differentiated/ 


14/16 


0.050 


0.669 


poorly differentiated 








Chemotherapy history/no 


25/5 


0.309 


0.61 1 


chemotherapy history 









Abbreviations: TNM, tumor-node-metastasis; KPS, Karnofsky performance status; 
LN, lymph node. 



Currently, a series of issues remain to be clarified. For 
example, it is suggested that a large number of IncRNAs 
detected in this study, ie, those other than AK124939, 
ASO3690, ENST00000491282, andXLOC_000575, are aber- 
rantly expressed in pulmonary adenocarcinoma tissues. Thus, it 
remains to be determined whether these IncRNAs are strongly 
correlated with survival in pulmonary adenocarcinoma and if 
they can be used as biomarkers for early clinical diagnosis. 

Clearly, there is much work to be done. RNA pull-down 
assays are indispensable both in vitro and in vivo. Moreover, 
describing the molecular mechanisms by which IncRNAs 
function in pulmonary adenocarcinoma is critical for elucida- 
tion of their role in cancer occurrence and development. Thus, 
IncRNAs have potential use as diagnostic and prognostic 
markers in pulmonary adenocarcinoma and in the develop- 
ment of IncRNA-mediated therapy. 

Conclusion 

We detected IncRNAs that are differentially expressed in pul- 
monary adenocarcinoma tissue compared with paired adja- 
cent noncancerous tissue from Chinese patients. AK1 24939 
was downregulated in pulmonary adenocarcinoma compared 
with nontumor tissues. Moreover, our results indicate that 
AK124939 levels in poorly differentiated adenocarcinoma 
tissue are lower than those in well to moderately differentiated 
adenocarcinoma tissue. In addition, gene coexpression net- 
works revealed a positive correlation between PDZD2 mRNA 
and AK 12493 9. Taken together, these IncRNA expression 
profiles in Chinese patients with pulmonary adenocarcinoma 
indicate that IncRNAs have potential use as diagnostic and 
prognostic markers in pulmonary adenocarcinoma and in the 
development of IncRNA-mediated therapy. 
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